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Abstract 

We propose and study a hydro dynamic model for pattern formation in mixtures of molecular 
motors and microtubules. The steady state patterns we obtain in different regimes of parameter 
space include arrangements of vortices and asters separately as well as aster-vortex mixtures and 
fully disordered states. Such stable steady states are observed in experiments in vitro. The sequence 
of patterns obtained in the experiments can be associated with smooth trajectories in a non- 
equilibrium phase diagram for our model. 



* E-mail : sumithra@imsc.res.in 

t E-mail : menon@imsc.res.in 

^ E-mail : suiiil@physics.iitm.ac.in 



1 



The mitotic spindle is a remarkable self-organized structure formed when a eukaryotic cell 
divides|l|. It consists of two separated inter-penetrating radial arrays (asters) of long, semi- 
fiexible polymeric filaments with polar character, called microtubules. Molecular motors 
such as kinesins and dyneins are protein molecules which "walk" unidirectionally along 
microtubules and exert forces on them, while consuming energy derived from the hydrolysis 
of ATP. Motor-microtubule interactions play a fundamental role in the formation of the 
spindle 



Interestingly, cell fragments containing bot 
to exhibit self-organized aster-like structures [2 



1 motors and microtubules have been found 
. Asters and vortices are seen in vitro, in 



experiments on mixtures of molecular motors, microtubules and ATP in a confined quasi- 
two-dimensional geometry Studying pattern formation in such simplified contexts may 
be useful in understanding the vastly more complex problem of spindle format ion 0, Q, 1^. 

These patterns form at large densities of motors and microtubules, rendering first- 
principles molecular simulations unviable. Such simulations require ab-initio modelling of 
interaction potentials; small errors in these potentials could be amplified at larger scales, 
thus changing the behavior fundamentally. Further, cellular pattern formation occurs un- 
der non-equilibrium conditions. These considerations motivate hydrodynamic approaches 
in which molecular scale information enters only in terms of coefficients in coarse-grained 
equations of motion for a few relevant fields P. Iffl]. 

An early attempt at describing such pattern formation is due to Lee and Kardar^^ and 
motivates the approach described here. Also, detailed Brownian dynamics simulations on 
a simplified model for motors, microtubules and their interactions have been performed by 



m 



12j |. However, these approaches 



Nedelec and collaborators, yielding valuable insights[ 
are unsatisfactory in some respects. The Lee-Kardar model fails to reproduce the complex 
patterns (arrays of asters, vortices and spirals) seen in experiments where the effects of 
boundaries are minimal. Further, it predicts that the unique steady state at high motor 
densities should be a single vortex whereas experiments obtain a lattice of asters in this 



limit 
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I3L Il4| . The simulations o: 



Nedelec et. al. are computationally intensive, requiring 
many parameters to be specified ja]. Several of these can be varied substantially without 
affecting the patterns formed - it is unclear which of them are crucial to pattern formation 
and which others play a secondary role. 

In this Letter, we propose a theory of pattern formation in mixtures of molecular motors 



2 



and microtubules. We motivate, using symmetry arguments, hydrodynamic equations of 
motion for a coarse-grained 2-dimensional vector field (T) representing the local orientation 
of microtubules as well as for local motor density fields rrif and rrih governing densities of 
"free" and "bound" motors. Motors which move on microtubules are called bound motors, 
whereas those which diffuse in the solvent are referred to as free motors. Bound and free 
motors interconvert at rates 7/-^6 and 7fe^/. Our equations describe the orientation of 
microtubules by complexes of bound motors, resulting in the formation of patterns at large 
scales. 

Before describing our model we briefly summarize our results: Our model generates 
virtually all the phases seen in the experiments, including a lattice of asters, a lattice of 



vortices /spirals 15l|. a mixture of asters and vortices as well as fully disordered phases fiol. 
We present a non-equilibrium "phase diagram" for our model; smooth trajectories in this 
phase diagram can be related to the sequence of states obtained experimentally as the motor 
density is increased. Such trajectories rationalize the difference in the sequence of patterns 
formed by different motor species - such as conventional kinesins and NCD's. We also 
predict the decay of bound and free motor densities in aster and vortex configurations. 

Our equations are derived in the following way: In the absence of interconversion terms 
changing a bound motor to a free motor, rrib obeys a continuity equation involving the 
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17|. The free motor field rrif 



current of motors transported along the microtubules 
obeys a diffusion equation with a diffusion constant D. These fields are coupled through 
interconversion terms 7/^?, (free — >• bound) and 7b^/ (bound free)j[]j|. We ignore spatial 
variations in the density of microtubules, concentrating on their orientational degrees of 
freedom. Our equation for T includes terms which reflect the alignment of microtubules 



by bound motors ll|. We also allow for a motor- independent stiffness against distortions. 
Finally, and most crucially for our purposes, symmetry permits a term of linear order in 
gradients of the bound motor density field to appear in the equation of motion for the 
microtubule field. Putting these ingredients together, we obtain 

dtuif = DV'^nif - 'if^bTrtf + Jb^frrih (1) 
dtiTih = -AV.{mbT) + 'jf^biTif - -fb^ffni, (2) 
dtT = T{a - (3T^) + -frubV^T + 7'Vm;,. VT + kV^T + SVrrib (3) 

The first term in the equation for describes the motion, with velocity A, of bound 



motors along microtubules. While 7/^6 should, in principle, depend on the local density 
of microtubule jl^. \v\ . we argue that such density variations can be neglected following 
coarse-graining to the mesoscopic length-scales of relevance to this problem and at the rela- 
tively large concentrations of microtubules in the experiments. We work in two dimensions 



throughout and scale length in units of D / AJ jS/a, time in units of /3D/{aA ), motor den- 



sities in units of D/7 and the tubule density in units of ya/(3. We choose 7' = e^y and vary 
the parameter e. 

The term in S can be thought of as arising from the variation of a free-energy-like term of 
the following kind: F ~ Jd'^x{{V.T) + Sf{mh)y. At non-zero S, such a term clearly favors 
a spontaneous divergence in the tubule configuration. When motor complexes bind to two 
parallel tubules, they pull these tubules inwards, generating a wedge-like configuration. Such 



12 1 which obtain a lattice of asters. This 



microscopic dynamics is incorporated in simulations 
term arises purely from bound motors and should, for small m^, have a strength linear in 
TUb i.e. finib) ~ nib. 

Three principal parameters enter our scaled equations: e, S and the total motor density 
m = rrif + nib. (We take 7/^b = 7b-^/ = 0.5 throughout and choose k to be 10% of m 
so that the motor- induced stiffness is dominant.) We solve the scaled system of equations 
numerically on an L x L square grid(l8|. The T equation is differenced through the Alternate 
Direction Implicit (ADI) operator splitting method in the Crank-Nicholson scheme while the 
equations of motion for the motor density are differenced through an Euler scheme. The 
grid spacing is 6x = 6y = 6 = 1.0, while the time step is 0.1. We evolve the equations for 
typically 2 x 10^ time steps, checking for convergence to steady state. The conservation of 
total number of motors is imposed by the boundary condition that no motor current flows 
into or out of the system. We choose "reflecting" boundary conditions for the T field, with 
all tubules at the boundary pointing along the inward normal jioj]. We work with systems 
of several sizes, ranging from L = 30 to L = 100, varying total motor densities in the range 
0.01 to 0.5 in dimensionless units [2^. 

Figures 1(a) - (d) depict four stable configurations obtained in different regimes of param- 
eter space for an L = 100 lattice. Fig. 1(a) shows a disordered arrangement of microtubules 
obtained at very low motor densities (m = 0.005) with e = 0.5 and S = 0. Figure 1(b) 
shows an aster-vortex mixture obtained at m = 0.01 at the same values of e and 5*. Note 
the presence both of well-formed asters and of vortices in the configurations. This figure is 
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to be contrasted to Fig. 1(c), obtained at m = 0.05, taking e = 5 and 5* = 0.001. Note 
the the absence of asters in this regime of parameter space. Finally, Fig 1(d), obtained with 
m = 0.5, e = 1 and S = 1, illustrates a lattice of asters phase, with asters being the only 
stable defects present. We can vary the sizes and numbers of asters obtained in configura- 
tions such as the one shown in Fig. 1(d), by changing S. A larger S yields a large number 



of small asters, while smaller values of S yield a smaller number of large asters 



17|. 



We also predict complex motor density profiles in aster configurations [2 1\ . Our results 
here differ both from those of Lee and Kardar as well as those of Nedelec, Surrey and Maggs 
(NSM)^3|- Assuming steady state aster configurations T = — f and solving the time- 
independent equations in steady state yields free motor densities which are a combination 
of confluent hypergeo metric Kummer functions j^]. The bound motor density field satisfies 
drmf{r) = —rriblr). The solution obeys the appropriate boundary conditions (no current at 
the boundary) and satisfies the normalization constraint on m. We obtain density profiles 
for free motors as 

where p = ^(1 , ), and the inverse of the "decay length" ^ is defined as = 

i'rb^f-y/i'rb-,f+i"/f-.b) 

2 









2p-l 



. The decay length ^ and the exponent p depend on 7/^6 and 

Provided ^ is large compared to the scale of the aster (we estimate ^ > AOfim for typical 
parameter values), the decay of motor densities is governed principally by the power-law 
nif ~ 1/r^. With = 7?,-^/ = 0.5, we obtain p = 1/3, although we can generate, as 7/^b 
and 76-*/ are varied, exponents which vary continuously in the interval [0:0.5]. (Power-law 
decays of mj are also obtained in the NSM model Ref. 



We relate our scaled parameters to typical experimental values in the following way: The 



tubule density, scaled in terms of ya/(3, is chosen to be unity. The diffusion constant D 
is about 20yum^/s and A ~ Ifim/s, defining basic units of length and time as 20/im and 
20 seconds respectively. A tubule density of 1 implies that over a coarse-graining length of 
400/im^, there are around 400 microtubules, a value close to that used in the simulations P|. 
The physical length of our simulation box is about 200/im, which is in the experimentally 
relevant regime. Our choice for '~ff~,b and 7b^/ corresponds to physical rates of 0.005s~^ to 
0.05s^^, slightly smaller than those in the simulations [8|]; using larger rates does not affect 



our conclusions here. 

Our results are summarized in the qualitative phase diagram of Fig. 2 which shows phases 
in the three-dimensional space spanned by e, S and m. For S* = 0, we obtain a disordered 
phase at low motor density, which becomes a lattice of vortices at somewhat higher motor 
densities. Large values of e (e > 1) yield well-formed vortex-like configurations while small 
e yields structures better described as aster-vortex mixture states. At intermediate values 
of e and m, spirals rather than vortices appear to dominate. At large m, with S = and 
large e, a single vortex is obtainedjll]. 

For non-zero but small S, these phases appear to continue out of the S = plane but 
are rapidly replaced by a lattice of asters phase for larger S. A cut of the phase diagram of 
Fig. 2 at finite S yields two "phases" - a disordered phase at small m and a lattice of asters 
at larger m. We can thus understand the sequence of patterns formed upon increasing m 
in mixtures of kinesin constructs with microtubules in terms of a trajectory which begins in 
the S = {or S sufficiently small) plane in the disordered phase and transits between the 
aster-vortex mixture and lattice of vortices phases (both of which lie in this plane) as m is 
increased. As m increases further and the effects of the S term become important, such a 
trajectory moves out towards non-zero S, encountering the lattice of asters phase. 

We have also examined the effects of changing the motor processivity, a quantity propor- 
tional to the ratio of 7/^6 to 76_>/. Smaller values of this ratio are appropriate to molecular 
motors such as NCD. At = 0.005, 7f,^/ = 0.05, we find that the disordered regime in 
the phase diagram of Fig. 2 expands, so that at equivalent values of m the disordered phase 
occupies much of the domain associated previously with the lattice of vortices. Whereas 
kinesins follow the sequence disordered - lattice of vortices - aster vortex mixture - lattice of 
asters as the m is increased, a mixture of microtubules with NCD motors bypasses the lat- 
tice of vortices phase altogether, transiting directly from the disordered phase to the lattice 
of asters phase in the experiments!^. In terms of Fig. 2, the expanded regime of disor- 
dered phase for the NCD motor suggests that phases such as the lattice of vortices and the 
aster-vortex mixture may be inaccessible at the motor densities at which the experiments 
are done, since the effects of the 5* term might be expected to dominate at large m. 

In conclusion, we have proposed a hydrodynamic theory capable of generating the se- 
quence of complex phases seen in experiments on pattern formation in mixtures of molecular 
motors and microtubules and constructed a non-equilibrium phase diagram which contains 
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these phases. The sequence of phases seen in the experiments on varying motor densities can 
be reproduced in terms of smooth trajectories in this phase diagram. Our resuhs are also 
predictive, since we can argue that smooth trajectories in parameter space cannot connect 
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disjoint phases directly. Further details will appear elsewhere 
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FIG. 1: Steady state configurations in our model at different parameter values (see text): (a) 
Disordered phase obtained at very low motor densities [m = 0.005 e = 0.5 and 5 = 0]; (b) Aster- 
vortex mixture obtained at [m = 0.01, e = 0.5 and 5" = 0]; (c) Lattice of vortices at [m = 0.05, 
e = 5 and S = 0.001]; (d) Lattice of asters obtained at [m = 0.5, e = 1 and 5 = 1] 
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FIG. 2: Schematic phase diagram in the e, m and S plane, illustrating the regimes of parameter 
space in which the phases of Fig. 1 are seen. 
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